Small Quarkonium states in an anisotropic QCD plasma 
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We determine the hard-loop resummed propagator in an anisotropic QCD plasma in gen- 
eral covariant gauges and define a potential between heavy quarks from the Fourier transform 
of its static limit. We find that the potential exhibits angular dependence and that binding 
of very small quarkonium states is stronger than in an isotropic plasma. 



I. INTRODUCTION 

Information on quarkonium spectral functions at high temperature has started to emerge from 
lattice-QCD simulations; we refer to ref. 1] for recent work and for links to earlier studies. This has 
motivated a number of attempts to understand the lattice measurements within non-relativistic 
potential models including finite temperature effects such as screening [3]. A detailed discussion 
of the properties of the heavy-quark potential in the deconfined phase of QCD is given in ref. Q , 
which also provides a comprehensive list of earlier work. 

In this paper, we consider the effects due to a local anisotropy of the plasma in momentum 
space on the heavy-quark potential. Such deviations from perfect isotropy are expected for a real 
plasma created in high-energy heavy-ion collisions, which undergoes expansion. We derive the 
HTL propagator of an anisotropic plasma for general covariant gauges, which allows us to define a 
non-relativistic potential via the Fourier transform of the propagator in the static limit. We also 
estimate the shift of the binding energy due to the anisotropy. 



II. HARD-THERMAL-LOOP SELF-ENERGY AND PROPAGATOR IN AN 

ANISOTROPIC PLASMA 

The retarded gauge-field self-energy in the hard-loop approximation is given by 

^' ^ J (27r)3 dk" p-v + iej - ^ ' 

Here, = (l,k/|k|) is a light-like vector describing the propagation of a plasma particle in 
space-time. The self-energy is symmetric, 11^'^ {p) = 11''^ (p), and transverse, p^n^'^(p) = 0. 

In a suitable tensor basis the components of 11'^'' can be determined explicitly. For anisotropic 
systems there are more independent projectors than for the standard equilibrium case [5|. To work 
in general covariant gauges, we use a four-tensor basis developed in ref. Q and the self-energy can 
now be written as Uf"" = aA''" + (33'^'' + -^Ct"" + SD'''' . 

In order to determine the four structure functions explicitly we need to specify the phase-space 
distribution function. We employ the following ansatz: 

/(p) = /iso(v/p2+C(p-n)2) , (2) 
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where n is the direction of anisotropy in momentum space and we choose n = (0, 0, 1) in this paper. 
The parameter ^ is used to determine the degree of anisotropy. Thus, /(p) is obtained from an 
isotropic distribution /iso(|p|) by removing particles with a large momentum component along n. 

Since the self-energy tensor is symmetric and transverse, not all of its components are inde- 
pendent. We can therefore restrict our considerations to the spatial part of H^'' and employ the 
contractions in ref. {§1 to determine the four structure functions. We do not list the rather cum- 
bersome explicit expressions for the four structure functions a, (5, 7, and 6 here since they have 
already been determined in ref. 

The propagator iA[^^ is diagonal in color and so color indices will be suppressed. In covariant 
gauge, its inverse is given by 
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where u) = p ■ m and m" is the heat-bath vector which equals to (1, 0, 0, 0) in the local rest frame. 
A is the gauge parameter. Upon inversion, the propagator is written as 
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-Ap^'p" , (4) 



A^'^{p'-a~j){u;'-(3)-S^p'-{n-pf] , (5) 
and the four vector rt^ = (0, n). For ^ = 0, we recover the isotropic propagator in covariant gauge. 



III. HEAVY QUARK POTENTIAL IN AN ANISOTROPIC PLASMA 

We determine the real part of the heavy-quark potential in the nonrelativistic limit, at leading 
order, from the Fourier transform of the static gluon propagator, 

Vir,0 = -g'C'F J |^e^P-A°°(a; = 0,p,C) (6) 

Here, Cf is the color factor and the ^-dependent masses m^, m|, and are given in ref. Q. 

One may check some limiting cases Q. The isotropic Debye potential is reproduced if ^ = 0. 
On the other hand, when r or ^ —^ 00, the potential then coincides with the vacuum Coulomb 
potential. Generally, the integral in ([7|) has to be performed numerically. The poles of the function 
are integrable. They are simple first-order poles which can be evaluated using a principal part 
prescription. In Fig. [1] we show the potential in the region f = rniD ~ 1 for various degrees of 
plasma anisotropy, where is the Debye mass. One observes that in general screening is reduced, 
i.e. that the potential at ^ > is deeper and closer to the vacuum potential than for an isotropic 
medium. This is partly caused by the lower density of the anisotropic plasma. However, the effect 
is not uniform in the polar angle. The angular dependence disappears more rapidly at small f , 
while at large f there is stronger binding for r parallel to the direction of anisotropy. Overall, 
one may therefore expect that quarkonium states whose wave- functions are sensitive to the regime 
f ~ 1 are bound more strongly in an anisotropic medium. 
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FIG. 1: Heavy-quark potential at leading order as a function of distance (f = rmo) for r parallel to the 
direction n of anisotropy and r perpendicular to n. 

Left: the potential divided by the Debye mass and by the coupling, V = V / {g^CFmo)- Right: potential 
relative to that in vacuum. 



IV. BINDING ENERGY FOR SMALL STATES 



Based on the above results, one can determine the binding energy of bound quarkonium states. 
In particular, we will concentrate on the shift of the binding energy due to the medium. An analytic 
estimate can be obtained from non-relativistic perturbation theory (to first order) for states with 
a Bohr radius times Debye mass small as compared to the anisotropy parameter ^. 

For weak anisotropy, ^ <^ 1, and distances f ^ 1, we expand the potential as 



V{y) w V^acir) + aniD - AiTa^m^, 
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where a = g^Cp/^T^- If mu/aniQ <C 1, where mg is the quark mass, the vacuum potential domi- 
nates and we can use Coulomb wavefunctions to calculate the expectation value of the perturbation 
given by the last two terms in the equation above. Finally, the binding energy for such small-size 
states can be expressed as 
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Note that E^ac < 0. The relative shift of the binding energy due to the presence of the (weakly 
anisotropic) medium is therefore given by 
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The first term corresponds to the reduced binding due to screening by the medium while the 
second term is the correction due to the non-zero anisotropy. The restriction to mjj/amQ <C ^ is 
necessary because terms of order 0{mD /amq) inside the brackets have been neglected. 

The above analysis applies also to excited states, provided that their size (in units of the Debye 
length) is small as compared to ^. The shift of the binding energy for the n*'' state is 
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However, it should be clear that for realistic cases (i.e., temperatures, quark masses and coupling 
constant), our above assumption that molamQ <C ^ is too extreme. For intermediate-size states 
and general ^, we must solve exactly the 3d Schrodinger equation with the anisotropic potential. 
This work is in progress. 

V. DISCUSSION AND OUTLOOK 

We have determined the HTL gluon propagator in an anisotropic (viscous) plasma in covariant 
gauge . Its Fourier transform at vanishing frequency defines a non-relativistic potential for static 
sources. We find that, generically, screening is weaker than in isotropic media and so the potential 
is closer to that in vacuum. Also, there is stronger binding of the quark pairs in the anisotropic 
system. Our results are applicable when the momentum of the exchanged gluon is on the order of 
the Debye mass to^ or higher, i.e. for distances on the order of Ad = l/mn or less. 

Following the discussion of the quark potential model in ref. at short distances, there is a 
string contribution to the potential which is not calculable perturbatively. However, at sufficiently 
high temperature, the perturbative contribution dominates over the linear confining potential at 
the length scale A^i. Roughly, this holds when T is larger than 2Tc. In this case, our result is 
directly relevant for quarkonium states with wavefunctions which are sensitive to the length scale 
Xd ■ Conversely, for those states whose length scale is larger, one should sum the medium-dependent 
contributions due to one-gluon exchange and due to the string [3] . 

At very short distances, the contribution from the string can also be neglected and one is 
dealing with perturbed Coulombic states. However, for charmonium or bottomonium, the string 
contribution is in fact important and it will be interesting to determine their wave functions and 
binding energies from a potential which is a combination of our anisotropic potential and the string 
contribution. 
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